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Abstract Stereoscopic white-light imaging of a large portion of the inner he- 
liosphere has been used to track interplanetary coronal mass ejections. At large 
elongations from the Sun, the white-light brightness depends on both the lo- 
cal electron density and the efficiency of the Thomson-scattering process. To 
quantify the effects of the Thomson-scattering geometry, we study an inter- 
planetary shock using forward magnetohydrodynamic simulation and synthetic 
white-light imaging. Identifiable as an inclined streak of enhanced brightness 
in a time-elongation map, the travelling shock can be readily imaged by an 
observer located within a wide range of longitudes in the ecliptic. Different 
parts of the shock front contribute to the imaged brightness pattern viewed 
by observers at different longitudes. Moreover, even for an observer located at a 
fixed longitude, a different part of the shock front will contribute to the imaged 
brightness at any given time. The observed brightness within each imaging pixel 
results from a weighted integral along its corresponding ray-path. It is possible 
to infer the longitudinal location of the shock from the brightness pattern in an 
optical sky map, based on the east-west asymmetry in its brightness and degree 
of polarization. Therefore, measurement of the interplanetary polarized bright- 
ness could significantly reduce the ambiguity in performing three-dimensional 
reconstruction of local electron density from white-light imaging. 

Keywords: white-light imaging, Thomson scattering, polarization brightness, 
interplanetary shock 
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1. Introduction 



Interplanetary space is permeated with the supersonic solar wind flow from the 
Sun. The ubiquitous solar wind not only carries large- amplitude Alfven waves 
(Bel cher and Davis, 197l| |Li and Li, 2007| ) but also is a transmission medium 
for interplanetary coronal mass ejections (CMEs) QDryer, 1994| ) . CMEs, firstly 
discovered by white-light coronagraphs in the 1970s, are large-scale expulsions 
of plasma and magnetic field from the solar atmosphere. A typical CME carries a 
magnetic flux of 10 23 Maxwells and a plasma mass of 10 16 grams ( |Forbes et al, 2006] ) . 
During interplanetary propagation, a CME may interact with the background 
structures in the ambient solar wind, such as the heliospheric current sheet (HCS) 
and corotating interaction regions (CIRs). Moreover, at solar maximum, multiple 
CMEs successively launched from the Sun are likely to interact with each other. 
The interaction is generally significant, nonlinear, and irreversible. The resultant 
complex travelling structures, sampled later at high temporal resolution by in 
situ instrumentation at 1 AU, are identified as magnetic cloud boundary layers 
QWei et al, 20031), CM E-CIR events ( |Dal Lago et al, 2006[), CME-shock events 
dLepping et al, 1997D , CME-CME events ( |Burlaga, Behannon, and Klein, 1987 



Dasso et al, 2009), and so on. Interplanetary CMEs have been widely recog- 



nized as the primary driver of interplanetary disturbances and large geomag- 
netic storms [e.g., |Burlaga, Behannon, and Klein| 119871 | Gosling et al.\ 119911 



Lepping et aT\ [T9971 [Dal Lago et at] l2006f 

Optical heliospheric imagers (His) enable continuous monitoring of the evo- 
lution of such solar eruptions as they propagate through interplanetary space. 
Heliospheric imaging fills the observation gap between near-Sun coronagraph 
imaging and in situ measurements. Interplanetary structures are viewed by 
means of Thomson-scattered sunlight. The scattered light retains the same spec- 
tral characteristics as the incident sunlight on the basis of Thomson-scattering 
theory ( Billings, 1966 Howard and Tappin, 2009] ) . The effects of multiple scat- 
tering are negligible, because the solar wind is optically thin. The angle be- 
tween the Sun and a target, such as a CME, as viewed from an observer, 
is termed elongation e. The brightness difference between the Sun and a tar- 
get (to be observed) at a large elongation is many orders of magnitude. More 
specifically, large CMEs at e = 45° have optical intensities that are of or- 
der 10~ 14 B Q ( |DeForest et al, 201 1[ ), where B Q , the solar brightness unit, is 
the intensity of the Sun as a power per unit area (of the photosphere) per 
unit solid angle. To map interplanetary structures at large elongations, the 
Thomson-scattered signals must be separated from many other sources of dif- 
fuse light: background light from the Sun, zodiacal light, and starlight (c.f., 
Figure 3 from Harrison, Davis, and Eylcsj 120051 Figures 5 from |Eyles et al 



2009] Figure 1 from lJackson et all I2010|) . The Thomson-scattering signals are 
faint and transient, whereas the other signals are intense but, fortunately, stable 
QLeinert and Pitz, 1989[ ). So by successfully subtracting the background bright- 
ness, an optical sky map taken by an HI easily reveals interplanetary transients 
irradiated by direct sunlight (e.g., fDcForost ~£ aUI201ip . IDeForest et 011120111 
proposed an image processing procedure consisting of five major steps: stationary 
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background removal, celestial background removal (including cross-image distor- 
tion measurement), residual F corona removal, moving feature filtration in the 
Fourier plane, and conversion back to focal plane coordinates. Using this proce- 
dure, the solar wind at 1 AU can be imaged with a sensitivity of a few x 10~ 17 B e , 
compared to a background signal of a few xlO _13 B Q . Instrument specifica- 
tions for an HI require careful design that incorporates the stray-light impacts 
of the HI bus, HI appendages, and other instruments onboard the spacecraft 
dHarrison, Davis, and Eyles, 2005"l|Eyles et al, 2009[ [Jackson et al, 2010D . His- 



toric milestones in terms of spaceborne white-light imaging were achieved by the 
zodiacal- light photometers ( |Leinert et al, 1981| ) onboard the Helios spacecraft, 
the Large Angle and Spectrometric COronagraph (LASCO) ( |Brueckner et al., 1995| )| 
onboard the Solar and Heliospheric Observatory (SOHO), the Solar Mass Ejec- 



tion Imager (SMET) (Eyles et al, 2003) onboard the Coriolis spacecraft, and the 



His ( |Howard et al, 2008[ |Harrison et al, 2008| |Eyles et al, 2009[ ) onboard the 



twin Solar- TErrestrial RElations Observatory (STEREO) ( |Kaiser et al, 2008D 
The STEREO mission is comprised of two spacecraft, with one leading the Earth 
in its orbit (STEREO A) and the other trailing behind the Earth (STEREO B). 
Both STEREO A and B separate from the Earth by 22.5° per year. The HI 
instrument on each STEREO spacecraft consists of two cameras, HI-1 and HI-2, 
whose optical axes lie in the ecliptic. Elongation coverage in the ecliptic is 4° - 
24° for HI-1 and 18.7° - 88.7° for HI-2; The field-of-view (FOV) is 20° x 20° for 
HI-1 and 70° x 70° for HI-2; The cadence is 40 minutes for HI-1 and 2 hours 
for HI-2 ( |Eyles et al, 2009[ |Harrison, Davies, and Rouillard, 2009] ). With the 
launch of STEREO, a CME can be imaged from its nascent stage in the corona 
all the way out to 1 AU and beyond, including those that are Earth-directed 

dHarrison et al, 2008||Davies et al, 2009||Liu~et al, 2010||Lugaz et al.,10To"l|DeForest et al, 20TT| ).| 

Moreover, images from the STEREO/ HI-2 have revealed detailed spatial struc- 
tures within CMEs, including leading-edge pileup, interior voids, filamentary 



structure, and rear cusps (DeForest et al, 2011). The leading-edge pileup of solar 



wind material is observed as a bright arc in optical imaging and revealed as a 
shock sheath from in situ sampling. With the launch of the STEREO mission, 
the interaction between CMEs in the inner heliosphere can be simultaneously 



imaged from multiple vantage points (Liu et al, 2012 1. His arc planned as part 



of the payload on a number of future missions, including the Solar Orbiter and 
the Solar Probe Plus. 

White-light images from His are more difficult to interpret than those from 
coronagraphs. CMEs imaged by the STEREO/ His are often faint and diffuse. By 
stacking a time series of running difference brightness along a fixed spatial slice, 
often corresponding to the ecliptic, a time-elongation map (commonly called a J- 
map) is produced. In such a J-map, the signatures of faint propagating transients 
are highlighted, and the transient can be tracked. However, the interpretation of 
the leading edge of the transient track in a J-map can be quite controversial. By 
fitting its time-elongation profile extracted from the J-map created from single- 
spacecraft observations dSheeley et al, 2008] |Rouillard et al, 2008D , the radial 
velocity and direction of a CME can be estimated, assuming that the CME 
moves radially at a constant velocity. Such a fit can be performed independently 
for the two STEREO spacecraft, if the CME is observed by both. Moreover, if 
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it is observed by both spacecraft, the CME's position (and hence its velocity) 
and its propagation direction can be inferred as a function of time through- 
out its observation by direct triangulation QLiu et ai, 2010[ ) or a "tangent to a 
sphere" method flLugaz et ai, 2010 ). The determination of CME kinematics, in 



particular the propagation direction, can be ambiguous. The assumptions that 
one is always observing the same part of a CME can lead to large errors in the 
estimated height of the CME's leading edge. The structure of a CME affects 
the derivation of its kinematic properties, particularly at large elongations. As 
a CME propagates to increasingly large elongations, the location of its lead- 
ing edge changes. As stated by |Howard and Tappm] 120091 the interpretation 
of observations of the leading edge of a CME can be fraught with difficul- 
ties, especially at large elongations. Scattered sunlight in the inner heliosphere 
is not only determined by the local electron density, but also by Thomson- 
scattering theory. The Thomson-scattering geometry is defined by a sphere, on 
which the Sun and the observer are located on opposite ends of the diameter. 
The classical theory of Thomson-scattering, by Billings) 119661 was revisited by 
|Howard and Tappin[ IMP [H oward and Tappin| 120091 concluded that: (1) The 
result of the scattering efficiency contribution is to somewhat de-emphasize the 
importance of the Thomson-scattering surface; (2) The Observed intensity is 
distributed out to large distances from the Thomson-scattering surface; (3) 
Such a spread of observed intensity is more significant at larger elongations. 
A CME can, therefore, be readily imaged at large elongations. Using SMEI 
observations, IHoward et ail 120071 demonstrated that limb CMEs could be de- 
tected out to at least 0.5 AU. Furthermore, |Tappin and How ard, 2009] showed 
that CMEs could be phenomcnologically modelled as a modified bubble or shell 
structure. However, inferring the three-dimensional location of a CME becomes 
more difficult with increasing elongation. In optical imaging, the observed in- 
tensity is traditionally assumed to come from a so-called plane-of-sky (POS). 
The POS is defined as the plane containing the Sun and being perpendicular 
to the line-of-sight (LOS) from an observer. While such a simple assumption 
is reasonable for coronagraph imaging, it is invalid for heliospheric imaging 
( Ho ward and Tappin, 2009] | Jackson et ai, 2010 ) . An inappropriate application 
of the POS assumption to heliospheric imaging would result in a significant 
underestimation of CME mass. A CME is inherently three-dimensional in na- 
ture, and imaging observations only provide two-dimensional information as a 
result of LOS integration through the three-dimensional structure. Effects of the 
projection and the Thomson-scattering are likely to affect brightness tracks in 
time-elongation maps. Uncertainties resulting from such effects can be assessed 
and disentangled quantitatively using a magnetohydrodynamic (MHD) mod- 
elling as a digital laboratory. For instance, CMEs on 24 - 25 January 2007 
were modelled, and synthetic optical images compared in detail with obser- 



vations from the STEREO/ His, by Lugaz et al. 120091 Moreover, by forward 



modelling of a travelling shock, |Xiong et al. 120111 demonstrated that remote 



sensing signatures in coordinated white light and interplanetary scintillation 
(IPS) observations can be theoretically simulated. Numerical models are useful 
in interpreting the brightness patterns of CMEs tracked in optical sky maps, and 
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hence aiding our understanding of hcliospheric processes such as the propagation, 
evolution, and possible interactions of CMEs. 

In this paper, we simulate a shock propagating through the inner heliosphere 
using an MHD model. Subsequently, we synthesize and investigate the corre- 
sponding white light images that would be observed from 1 AU. We present 
the numerical MHD model and the optical Thomson-scattering formulation in 
Section [2] In Section [3j we describe tracking of the shock front from different 
perspectives, and present profiles of the scattered sunlight intensity along various 
LOSs. Subsequently, we analyze the brightness patterns observed in synthetic 
time-elongation maps and longitude-elongation maps (Section 0]). Finally, we 
discuss remote sensing of much more complex interplanetary phenomena in both 
white light and radio in Section [SJ 



2. Method 



A numerical model of the heliosphere QXiong et at, 2006a[ ) is used to associate 
interplanetary dynamics with corresponding observed signatures. Such a self- 
consistent link is summarized in the flow chart for synthetic remote-sensing 
in white light and radio shown in Figure 8 from Xiong et al.\ 120111 Here, we 
investigate only white-light imaging for a travelling fast shock. The procedure 
consists of two steps, described in Sections 12.11 and 12.21 



2.1. Numerical MHD Model 



Large-scale dynamics of the inner heliosphere can be physically described by 
MHD processes. Using a sophisticated numerical MHD model ( |Xiong et ai, 2006a[ ),[ 
a blast shock wave travelling through interplanetary space is numerically simu- 
lated in this paper. We describe our model in Table [TJ establish the background 
solar wind in Table [3J and prescribe the shock injection in Table [3] The simulated 
electron density is used to generate synthetic white-light images. 



2.2. Thomson-Scattering Calculations 



White-light imaging is performed on the basis of the well-established Thomson- 
scattering theory. Interplanetary free electrons scatter white-light photons from 
the photosphere, so the inner heliosphere is visible as a brightness sky map. 
In this paper, we adopt the formulation of Thomson-scattering theory in in- 
terplanetary imaging given by |Howard and Tappin 120091 This is demonstrated 
schematically in Figure[TJ At a scattering point Q, the solar surface is only visible 
within a cone where ASQT = Cl and ZSTQ = 90°. The Sun looks like a limb- 
darkened luminous disk. All straight rays of direct sunlight intercepted at point Q 
are scattered into a full solid angle of 4 it, with part of the sunlight being scattered 
toward an observer at point O. The observer O, scattering point Q, and Sun cen- 
ter S all lie in an emergent plane. At the observation site O, received sunlight is 
direct and unpolarizcd along the LOS SO, and scattered and polarized along the 
LOS QO. The luminous magnitude of the Sun (LOS SO) is dramatically higher 
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than that of interplanetary space (LOS QO) ( |Harrison, Davis, and Eyles, 2 005; 
|Eyles et al, 2009| |Jackson et al., 2010| |DeForest et al, 2011[ ). A series of baf- 
fles onboard the STEREO/ His (Eyles et al., 20091 attenuates the intense solar 
radiation to such a degree that the much fainter interplanetary brightness is 
revealed. Scattered by a single electron at point Q, the optical intensity G, 
measured at point O as a power per unit solid angle, can be expressed in units 
of B . The intensity G is polarized along the LOS QO, and can be expressed in 
terms of two orthogonal components Gr and Gt- The radial component Gr is 
in the previously defined emergent plane (the OQS'-plane) and is perpendicular 
to the LOS QO. The tangent component Gt is perpendicular to the OQS- 
planc. Mathematical expressions for the luminous intensity G and its constituent 
components (Gr and Gt) are formulated in terms of van de Hulst coefficients 
A, B, G, and D ( |Howard and Tappin, 2009) ): 



Gp — Gt — Gr 
G = Gt + Gr = 2 Gt — Gp 
Gt = 



G P = — sin x [(1 



£[(l-u)C + tt£>] 

u)A + uB] 



(1) 



Here Gp is an auxiliary parameter. a c , defined as the differential cross section for 
perpendicular scattering, is aconstant of 7.95xl0~ 30 m 2 Sr _1 (Howard and Tappin, 2009). 
u, defined as a limb-darkened coefficient, is a function of wavelength. A value 
for u of 0.63, for a wavelength of 5500 A is adopted in this paper, z is the 
distance from the scattering point to the observer. The coefficients A, B, C, and 
D ( |Billings, 1966[ ) are given as: 



A = cos SI sin 2 SI 
1 



B 



I - 3 sirT O 



cos si 

sin SI 



(1 + 3 sin 2 SI) In 



f + sin SI 
cos SI 



4 cos 3 n 
3 ~ °° 3 — 

cos 2 n 



(2) 



D = 



5 + sin SI — 



sin SI 



(5 -sin 2 SI) In 



f + sin SI 
cos SI 



As the scattering point Q increases in heliocentric distance r, ZSQT = SI 
becomes smaller and the Sun tends toward a point source. Accordingly, the 
coefficients A, B, C, and D fall off as a function of r~ 2 , and the differences given 
by \A — C\ and \B — D\ reduce. As a result, the incident sunlight becomes fainter 
and more collimated (c.f., Figure 4 from |Howard and Tappin| I2009p . 

An imaging detector collects photons that fall within its FOV, and the bright- 
ness within a pixel will comprise contributions from all scattering sites along its 



LOS. Howard and Tappin 2009] pointed out, however, that " When interpreting 
real observations, it is crucial to note that although the scattering is presented as 
a line-of-sight integral this is not strictly true. It is in fact an integral though the 
cone of the instrument's point spread function." . The configuration for the LOS 
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integral calculation is demonstrated schematically in Figure [TJd. The scattering 
plane is at a distance z from the detector. The detector itself has an area of 
<5A, and its FOV at distance z has an area given by 5ui = dx dy. The solid 
angle subtended by the detector at the scattering area 8u> at z is given by 
du> = Suj/z 2 = dxdy/z 2 . The brightness contribution from distance z to z + dz 
along the LOS, where the electron number density is given by n, can be expressed 
as follows: 



( * \ 
If we let 



dz = JT n Bq 

5uj 



( G\ 

G P 

Gt 
\G R J 



dx dy dz 6 A = n z 2 



( G \ 

Gp 

Gt 
\G R J 



Bp, du>SA = 1. 



B Q du> 6 A dz . (3) 



(4) 



Equation ([3]) can be further simplified to 

/ i \ 

ip 



( G \ 

Gp 

Gt 
\G R J 



(5) 



Within each pixel of the detector, the optical brightness that is recorded is the 
LOS integral of Equation ((SJ: 









Ip 




ip 


It 


Jo 


i T 


\Ir) 







dz 



Gp 
Gt 
\GpJ 



dz . 



(6) 



Hence it can be seen that the electron density n and Thomson-scattering ge- 
ometry factors (z 2 G, z 2 Gr,, z 2 Gt) jointly determine the observed brightness. In 
terms of the electric field oscillations 7r and It, polarization p is defined as: 

— Ir Ip 



P = 



h + In 



(7) 



Sunlight is scattered backward for \ < 90° (Figure QJ) , perpendicular for 
X = 90° (Figure QJl), and forward for \ > 90° (Figure QJ). Perpendicular scat- 
tering toward a fixed observer O comes from the so-called Thomson-scattering 
sphere, which is centered between the Sun and the observer (Figure QJl). On 
the Thomson-scattering surface itself, Thomson-scattering is actually minimized 
due to perpendicular scattering. However, both incident sunlight intensity and 
local electron density are maximized on the Thomson-scattering surface, that 
being where any LOS is closest to the Sun. These three competing effects de- 
termine the scattered sunlight intensity. Although maximized on the Thomson- 
scattering surface itself, the scattered intensity is distributed with distance away 
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from the surface ( |Howard and Tappin, 2009] ). Owing to the weakest scatter- 
ing efficiency being on the Thomson-scattering sphere, the importance of the 
Thomson-scattering sphere is de-emphasized, particulary at large elongations. 
On the Thomson-scattering sphere, the radial component Jr. of scattering is 
the smallest and tangent component It is the largest. Thus, the polarization 
of scattered sunlight is largest on the Thomson-scattering sphere, according to 
Equation ([7]). If it can be measured at large elongations, the polarization of 
scattered sunlight can provide an important clue to identifying corresponding 
major scattering location. 

3. Line-of-Sight Profiles of Scattered Sunlight 

The position of an interplanetary shock can be identified from a white-light 
imaging map as a locally enhanced brightness pattern. Sunlight scattered from 
a travelling shock is more intense than that scattered by the ambient solar 
wind, as the electrons are significantly compressed just downstream of the shock 
front. For the purposes of this paper, we consider an interplanetary shock viewed 
simultaneously by three observers that are located at fixed longitudes of ip = 0° , 
—45°, —90°. The observers and their corresponding Thomson-scattering spheres 
are depicted in Figure [2J The observed brightness pattern looks different when 
viewed from the three different sites, and its leading edge corresponds to different 
sections of the shock. As shown in Figure |31 the nose, left flank, and right flank 
of the shock are penetrated by LOSs LI, L4, and L7, respectively. Ray-paths 
LI, L4, and L7 are directed toward different observers, at longitudes of ip = 0°, 
—45°, —90°, respectively. Each of the ray-paths LI, L4, and L7 corresponds 
to the largest elongation at which enhanced brightness is observed in a full 
sky map. Even for a fixed observer, the leading edge of the brightness pattern 
actually corresponds to a different section of the shock at a given time. This can 
be demonstrated by considering an observer-directed shock travelling along the 
longitude <p = 0° (Figure [3]). The travelling shock front is intersected by LOS 
LI at 12 hours, L2 at 25 hours, and L6 at 37 hours, respectively. The shock 
flank contributes to the brightness pattern observed at 12 hours, and the shock 
nose at 37 hours. Such a transition from the flank to the nose is smooth in 
time. With respect to the observer, plasma packets along a single LOS generally 
have different spatial speeds (where a spatial speed is expressed in an observer- 
centered polar coordinate system). Different spatial speeds generally correspond 
to different observer-centered angular speeds. So, plasma packets temporarily 
aligned along one LOS will spread over adjoining LOSs later (c.f., Figure 6 from 
|Xiong et al.\ l20lTj) . It should be remembered that when one visualizes optical 
signatures of an interplanetary CME in a time-elongation map, he/she is not 
always observing the same part of the structure. Interpretations of brightness 
patterns at large elongations are usually ambiguous. 

Sunlight scattered along an LOS is collected in the corresponding pixel and 
recorded in terms of optical brightness. The brightness in the pixel results from 
an integral along the LOS, as expressed in Equation (jB]). Thus, the LOS dis- 
tributions of localized intensities (i, ir, zt) and Thomson-scattering geometry 
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factors (z 2 G, z 2 Gr, z 2 Gt) will determine which part of the LOS contributes 
most significantly to the total intensities observed (I, 7r, It)- Such analysis of 
an interplanetary shock is performed for LOSs L1-L7 in Figure[3l and the salient 
parameters are presented as a function of modified scattering angle x* = 90° — X 
in Figures 0] and [5J Note that Figures @] and [3J present normalized values of 
the intensities, Thomson-scattering factors, and electron density. The normal- 
ization factors are given in Table El The Thomson-scattering factors depend 
on the relative geometry between the Sun, a scattering site, and the receiving 
observer. More specifically, the elongation e and scattering angle x determine the 
Thomson-scattering factors. The dependence of Thomson-scattering factors on 
e and x is explored by |Howard and Tappin 120091 For a fixed e, the distribution 



of the Thomson-scattering factors is symmetric around x* = 0° (Figures Hfc, HJ, 
E£, 0, andJSJa). x* = 0° corresponds to perpendicular scattering, i.e., x = 90°, 
and the scattering site lies on the Thomson-scattering sphere itself (Figure QJL). 
For 90° < |e| < 180°, an observer views the hemisphere opposite the Sun, and all 
optical paths toward the observer are backward with x* > 0° . One such example 
is LOS L6 in FigureUJ The profiles of the Thomson-scattering factors are convex 
around x* = 0° f° r z 2 Gt and concave around x* = 0° f° r z2 Gr (Figures Ht, HI, 
[5b, [5J, and [5b). At large elongations, z 2 Gr is negligible at x* = 0°. Moreover, 
at x* =0°, the scattered sunlight is nearly linearly polarized, and its polariza- 
tion p is almost 1. As elongation e increases, the profile of z 2 G — x* flattens. 
Along the nearly anti-sunward LOS L6, where e = 155°, the scattered sunlight 
is almost unpolarizcd with Gr f» Gt and p f» (Figures 0Ll and HJjv). The 
Thomson-scattering factors remain fairly constant within that region delimited 
by X*l < 30°. In Figure [31 an incident shock compresses plasma at its front, 
leaving a low density void behind. The distributions along LOS L3, which cuts 
across the density void and the shock flanks (Figures [3b and [HJi), are presented 
in Figures [5ji-f. Figure [5ji reveals density spikes within —49° < x* < —40° 
and 25° < x* < 35°, and a density void within -40° < x* < 25°. The 
locally scattered intensity is jointly determined by the electron density and the 
Thomson-scattering factors, with i = nz 2 G, in = n z 2 Gr, and it = nz 2 GT- 
The LOS profile of scattered intensity i (Figure [5Ji) is similar to that of electron 
density n (Figure [5ji) ; this is not the case for the background signals. As it is 
7, the integral of i along the LOS that is recorded by a corresponding pixel, 
the two density spikes in Figure [SH would be essentially undetectable, based 
on measurements along LOS L3. Other LOSs in Figure [3] are, however, tangent 
to the shock front. For an observer situated at a longitude ip = 0°, the trav- 
elling shock can be continuously tracked as a moving brightness pattern. At 
12 hours, when the leading edge of the pattern is observed at an elongation of 
e = 17.5°, the shock front is bounded between —49° < x* < —31° along this 
LOS, LI (Figure |4ji). Simultaneously, a spike-like density enhancement exists 
at x* = 15°; This spike, which is the signature of the heliospheric plasma sheet 
(HPS), has shifted from its initial location of x* = 8°. The HPS, initially located 
at if = 90° and with a longitudinal width of only 3°, co-rotates with the Sun. 
The density of HPS reduces with increasing heliocentric distances, so its optical 
brightness is negligible at large elongations. Thus the travelling shock contributes 
most significantly to the observed transient brightness enhancement. As time 
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elapses, the observer-directed shock is viewed at increasingly larger elongations, 
and moves nearer to the corresponding Thomson-scattering sphere. Accordingly, 
region inan-x* plot that corresponds to the shock compression gradually shifts 
toward x* = 0° (Figures HJi,|l], and Hp). Once the observer-directed shock front 
has passed over the observer at 1 AU (Figures [3^ and[3j), it scatters sunlight 
back toward the observer. For such backward scattering, the Thomson-scattering 
factors decrease monotonically with an increase in LOS depth, z (Figures 0p 
and Hli). For an observer viewing along elongations e > 90°, electrons in the 
vicinity of the observer contribute to the majority of observed total brightness, 
/. Hence optical imaging at elongations e > 90° can provide measurements of 
electron number density in a similar manner to in situ observations made by 
a spaceborne electron detector. For an observer situated at <p ^ 0°, the shock 
is propagating off the Sun-observer line. However, the shock front is very likely 
to cross the observer's corresponding Thomson-scattering surface. In this case, 
the most favorable configuration for white-light imaging is when the nose of the 
shock is on the Thomson-scattering surface, as can be seen from LOS L4 in 
Figures and [Sg-1. The wide flank of the shock wave is readily imaged, even 
if the shock is propagating perpendicularly to the Sun-observer line. Imaged 
from tp = —90° (Figures [5^ and[5f), the nose and left flank of the shock are 
too far away to be detected whereas the right flank of the shock can be detected 
along LOS L7 (FigurcsjSJii-r). At large elongations, both the Thomson-scattering 
sphere and the CME volume must, therefore, be considered when interpreting 
optical imaging of the inner heliosphere. 

4. Optical Brightness Patterns 

Moving brightness patterns in optical sky maps provide observational evidence 
of interplanetary transient disturbances such as travelling CMEs. CMEs often 
travel faster than the ambient solar wind, resulting in compression. Thus CMEs 
are often associated with dense plasma, such that sunlight scattered within the 
CME volume is intense. The wide FOVs of the new generation of HI instru- 
ments encompass a wide range of elongations. Transient brightness dramatically 
reduces with an increasing elongation e. However, the observed brightness over 
the entire FOV can be rescaled to be within the same order of magnitude, using 
the normalization factors shown in Figure [6j In the generation of Figures [6^-d, 
the electron number density n is assumed to vary with heliocentric distance r 
according to n oc r~ 2 . Intensity profiles (z, ir, it) are then calculated using the 
same method as for Figures H] and [5] Assuming that n oc r~ 2 , and using the 
known values of the Thomson-scattering geometry factors (z 2 G, z 2 Gr, z 2 Gt), 
integrated intensities (/*, I T ) within 0° < e < 180° are derived; These are 
presented in Figures Eb-g. /*, and I T , which are linearly scaled to be 1 at 
e = 6.7°, are used to normalize the intensity I, Jr., and It shown in Figures [7] 
and[8l CMEs cause significant derivations in heliospheric electron density n away 
from the equilibrium situation dictated by n oc r~ 2 . These normalized intensities, 
I/I*, Ir/Ijl, and It/I t , can be used to identify transient features in brightness 
sky maps. In this paper, we present such normalized intensities, for the case of 
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a travelling shock (Figure [3]), as longitude-elongation {tp — e) maps (Figure [7]) 
and time-elongation (t — e) maps (Figure E]). In Figures [7Jt-i and |5h,-i, regions 
with I/I* > 3.68 x 10~ 15 , I R /I R > 0.68 x 1CT 15 , and I T /I% > 2.93 x 10~ 15 
arc delimited by dotted lines, and are identified as corresponding to the shock 
front. The mapped brightness pattern significantly depends on the longitudinal 
position of the observer (Figure [7]) . An observer at a location within the range 
150° < \tp\ < 180° is nearly opposite to the direction of shock propagation, so 
optical signals of the shock arc too faint to be detected. An observer at tp = 0° 
cannot detect the shock until t = 12 hours. The patterns of brightness within 
the longitude-elongation maps are symmetrical with respect to elongation for 
ip = 0° and unsymmetrical elsewhere (Figure [7J. The propagating shock, as 
observed from a fixed longitude, is manifested as an inclined streak in a time- 
elongation map (Figure [5]). The slope of this feature is steepest for <p = 0° 
(Figure [5]), as its arrival time at e = 100° is 35 hours at ip = 0° compared with 
50 hours at tp = —45° and 90 hours at tp = —90°. The most intense optical 
signatures of the shock are those observed from tp — ±60° at a time of 12 
hours (Figure UJi) and from tp = ±30° at 25 hours (Figure [7}d). From these 
vantage points, the brightness patterns corresponding to the shock cover the 
widest elongation extent and exhibit the largest relative brightness enhancement 
due to the close proximity of the shock nose to the Thomson-scattering surface 
at these times (see, for example, L4 in Figure Efc). When the shock front crosses 
1 AU, the optical brightness is enhanced over the entire range of elongations 
90° < |e| < 180°, simultaneously. Such a crossing occurs at around t = 37 hours 
for an observer at tp = 0° (Figures [Jj: and [8^,), and near t — 53 hours for an 
observer at tp = —45° (Figure [SJd). The brightness enhancement, which results 
from backward scattering of sunlight, lasts for between 15 and 20 hours (Figures 
[5^-b). These figures demonstrate how easily an interplanetary shock, with its 
wide front, can be optically imaged. 

Polarization measurements of scattered sunlight are very useful in locating 
the three-dimensional position of a volume of dense plasma. Viewed along any 
elongation, the polarized intensities (/, /r. It) are an LOS integral of the 
local electron density, weighted by Thomson-scattering geometry factors. The 
Thomson-scattering factors arc different for the different polarized components, 
i.e., z 2 Gr for a radial component ir and z 2 Gt for a tangent component %t in 
Equation ©. Near the Thomson-scattering surface (where x* ~ 0°), z 2 Gt is 
much larger than z 2 Gr (Figures Hb, HJ, Hp, [St, 0, and[5b), and the locally scat- 
tered sunlight is nearly linearly polarized. When the heliospheric electron density 
n is in equilibrium, the distribution of n depends only on heliocentric distance 
r, according to the expression n oc r~ 2 , and the integrated scattered sunlight 
along elongation e = 70° has the largest polarization, with a value of almost 0.8 
(Figure Hh) ■ An observer viewing sunlight scattered by a dense parcel of plasma 
will detect an increase in /, an increase in It, a decrease in Jr, and an increase 
in polarization p as the parcel of plasma approaches the Thomson-scattering sur- 
face. Polarization measurements for the shock studied in this paper are presented 
as longitude-elongation maps in Figure [7] and as time-elongation maps in Figure 
HI LOSs L2, L4, L5, and L7 (in Figure EJ) have polarizations of 0.68, 0.8, 0.81, 
0.7, respectively. Obviously, when the shock nose lies on the Thomson-scattering 
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surface itself, the polarization of scattered sunlight is significantly enhanced. In 
the longitude-elongation (ip — e) maps shown in Figures 0-1, two localized and 
symmetrical patches of enhanced polarization are evident. At an elapsed time of 
25 hours, one of these patches lies within the region defined by 40° < e < 55° 
and —75° < ip < —10° (Figure [TJc) . Inferring the three-dimensional location 
of an interplanetary CME at large elongations using a total intensity measure- 
ment is fraught with uncertainties ( |Howard and Tappin, 2009[ ), which could be 
significantly reduced through the use of additional polarization measurements. 
Locating a CME using measurements of polarization would require the theo- 
retical calculation of the Thomson-scattering geometry factors within the CME 
volume. However, once the inherent effects of Thomson-scattering are removed, 
the total optical intensity can be used to reliably estimate the mass of the 
CME. A full sky imager observing polarized light would not only be capable 
of monitoring the inner heliosphere, but would also enable the three dimensional 
location of interplanetary disturbances to be established. 

The characteristic signature of a large-scale interplanetary transient in a 
brightness sky map is a bright arc, followed by a dark void. The bright arc 
corresponds to an interplanetary shock. Such a shock can be further classified 
either as a blast shock or a piston-driven shock. A blast shock is generated in 
response to an impulsive release of energy in the low corona such as a solar 
flare; a piston-driven shock is usually formed ahead of a fast CME. Coronagraph 
observations of loop-like transients are inconsistent with theoretical models of a 
blast shock wave (Sime, MacQucen, and Hundhauscn, 1984). Here, wc consider 
only the white-light signatures of the interplanetary shock; the exact nature of 
the shock is not of relevance. In practice, features in a brightness sky map need 
to be detectable in absolute intensity and exhibit a sharp gradient along their 
boundary. Satisfying the criteria of detectability indicates that an interplanetary 
disturbance should be relatively close to the Thomson-scattering sphere, i.e., 
usually within \x*\ < 30° (Figures [TJ;-e). However, the shock studied in this 
paper is very strong, described by initial parameters given in Table [3] If the 
shock was much weaker, the angular extent of its front would be less. In the 
case of a weak shock, for an observer at ip = 0°, the shock front would lie 
well inside the Thomson-scattering sphere and would be completely invisible 
before arriving at 1 AU. A bright arc observed by an optical imaging device may 
correspond to the shock forerunning a CME. During their propagation, CMEs 
continuously change their position relative to the Thomson-scattering surface 
of a corresponding observer. Therefore, the white light signature of a travelling 
CME will evolve continuously. 



5. Discussions and Summary 

In this paper, we present the results of a forward modelling study of an incident 
shock in the inner heliosphere. Through combining a numerical MHD model 
( Xiong et ai, 2006a ) and Thomson-scattering theory ( |Howard and Tappin, 2009] ) .[ 
we are able to investigate the causal link between interplanetary dynamics 
and observable signatures. By synthesizing white-light observations, we confirm 
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the significant role played by the Thomson-scattering geometry in determining 
the brightness of sky maps at large elongations. Furthermore, we suggest the 
use of interplanetary polarization measurements to help locate interplanetary 
CMEs; No such interplanetary observations are currently made. In contrast, the 
STEREO/ COR and LASCO coronagraphs do have polarizers, and have been 
used successfully to investigate CME orientations, close to the Sun, using a po- 
larimetric reconstruction technique (e.g.. [Moran et 01112010]) . As demonstrated 
by the proof-of-concept study presented here, however, interplanetary polariza- 
tion measurements made possible by future technological advances would be of 
significant benefit in the study of interplanetary transients such as CMEs. 

Optical brightness maps contain imprints of various large-scale interplanetary 
processes. These processes are usually much more complex than a simple case 
of a travelling shock, as studied in this paper. The use of numerical models 
is crucial in uncovering the physics that underlies white-light observations, such 
that optical imaging technology can be exploited to its fullest. To our knowledge, 
We list some examples below: 

1. Interplanetary coupling of multiple CMEs/shocks. 

Successive CMEs can interact with each other during their propagation, and 
form compound cjecta within 1 AU. A high-speed CME can drive a fast shock 
ahead of itself, which is wide in angular extent than the CME body itself. 
So CME-shock interaction is much more likely than CME-CME coupling. 
Both CME-shock interaction ( |Xiong et al., 2006a] |Xiong et al, 2006b] ) and 



CME-CME coupling QXiong et al, 2007[ |Xiong et al., 2009[ ) were theoreti- 



cally simulated in terms of interplanetary dynamics and ensuing geoeffec- 
tiveness; These theoretical modelling results are observationally confirmed by 
a recent work based on STEREO/ HI imaging of the interaction of multiple 
interplanetary CMEs ( |Liu et al., 2012[ ). The most compelling evidence for 
CME-CME interaction in an optical sky map is the formation of an intensely 
bright arc as the CME-driven shocks completely merge. This bright arc cor- 
responds to the merged shock front, a region in which the plasma is greatly 
compressed. Photosphcric data can be used as a driver for numerical MHD 
models of the inner heliosphcre, to enable realistic simulations of multiple 
CMEs to be performed. Following the same method as used in this paper, 
optical sky maps can be synthesized and then compared with observed sky 
maps. Such a comparison helps validate the numerical MHD modelling. The 
use of validated numerical modelling enables individual CMEs and interaction 
regions to be identified within the complex brightness patterns that can be 
present in optical sky maps. 

Interplanetary coupling of a CIR with CMEs/blobs. 

CIRs are periodic spiral-like structures formed in the inner heliosphere as a 
result of compression at the interface between fast and slow streams. CIRs 
form where the fast solar wind is emitted from the rotating Sun, behind the 
slow solar wind, along the same solar radial. Like CMEs, although smaller 
in scale, plasma blobs contribute to the more transient nature of the inner 
heliosphere. Coronagraph observations show that plasma blobs are intermit- 
tently released from the cusps of coronal helmet streamers, and propagate 
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at the typical slow solar wind speed ( jWang et al., 1998] ). A CIR can sweep 
up slow blobs in front, or can impede the propagation of a following fast 
CME; Such blobs and/or CMEs are entrained by the CIR. Plasma blobs 
that become entrained within the CIR undergo strong compression such that 
they are highly visible in optical sky maps. Blobs entrained at the stream 
interface can be identified as recurring patterns of inclined streaks in time- 
clongation maps. In the case of CIR-CME interaction, a CIR could be warped 
by the entrained CME. CIR-CME interaction is generally much more dis- 
ruptive than interaction between CIRs and blobs. Optical observations of 
blobs that have become entrained at the stream interface have been pre- 
sented by |Sheeley et al\ 120051 and IRouillard efoO 120051 CIRs are inferred 
to have the spatial morphology of a garden-hose density spiral, based on 
their signatures in time-elongation maps. Interplanetary interaction between 
CIRs and CMEs/blobs can result in complex brightness patterns in optical 
sky maps. The observed brightness is dictated by multiple factors such as 
the Thomson-scattering geometry, the alignment of the spiral CIR along the 
LOSs, and local compression within interacting sites. Numerical models can 
aid our understanding of these CIR-rclatcd phenomena. 
3. Coordinated observations from multiple optical and radio telescopes. 

Simultaneous heliosphcric imaging from multiple vantage points was suc- 
cessfully realized with the launch of the twin STEREO spacecraft. With 
stereoscopic imaging, multiple viewing ray-paths from one observer intersect 
with those from the other observer. Thus the three-dimensional distribu- 
tion of electrons in the inner heliosphere can be reconstructed using a time- 
dependent tomography algorithm fc.f.. Uackson et aOl2003llBisi et aZ.l I2010|) . 
At times, several IPS ray-paths lie within the imaging FOV such that CMEs 
can be simultaneously observed in IPS and white light ( |Dorrian et al., 2008| 
|Manoharan, 2010] ). The IPS technique can be used to estimate the location 
and speed of micro-scale electron density irregularities ( Hewish, Scott, and Willis, 1964| | 



Coles and Harmon, 1989 ). If such irregularities exist within a CME and some 
assumptions are made about the CME kinematics, the IPS signals can be used 
to predict the appearance of the CME in later optical sky maps ( |Xiong et al., 201lj ).| 
It should be noted that newly constructed low-frequency radio interferometers 
such as the LOw Frequency ARray (LOFAR) ( |de Vos, Gunst, and Nijboer, 2009"1 )| 



Lonsdale et al., 2009 1 are prov 



and the Murchison Widefield Array (MWA) 
ing a major milestone in IPS technology. Coordinated remote sensing ob- 
servations in the optical and radio regimes enable the inner heliosphere to 
be mapped in fine detail; Numerical models can guide such endeavors by 
suggesting synthesized results beforehand (e.g., |Xiong et aZ.(l201ip . 

In closing, white-light imaging is a mainstream technology for remotely sens- 
ing the inner heliosphere. However, owing to inherent geometry effects of Thomson-^ 
scattering, deriving the kinematic properties of interplanetary transients from 
optical sky maps is ambiguous at large elongations. The ambiguities arising 
due to the Thomson-scattering geometry can be rigorously constrained, if op- 
tical imaging is complemented by other observational techniques, such as radio 
imaging and in situ sampling, and/or by numerical modelling. More theoretical 
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modelling will be done as a natural extension to the preliminary results presented 
in this paper. 
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Table 1. Description of our numerical heliosphere model (Xiong et al., 2006a) 



Mathematical description 


MHD equation set 


scheme 


Total Variation Diminishing (TVD) + 




Weighted Essentially Non-Oscillation (WENO) 


type 


2.5 dimensional 


domain 


the ecliptic 


coordinate 


radial distance r, longitude ip 


boundary 


25 < r < 230 solar radii, -180° < <p < 180° 


mesh size 


Ar = 0.5 solar radii, Aip = 0.5° 



Table 2. Prescription of a background solar wind at the inner 
boundary of 25 solar radii. 



type 


slow solar wind 


structures 


spiral interplanetary magnetic field, 




a hcliosphcric plasma sheet astride 




a heliospheric current sheet (HCS) 


initial longitude of the HCS 


V = ±90° 


species 


proton p, electron e 


number density 


n = n p = n c = 550 cm -3 


radial speed 


v r = 375 km s _1 


magnetic field strength 


B = 400 nT 


temperature 


T p = T c = 9.6 x 10 5 K 


plasma beta 


= 0.23 


alignment condition 


V || B 
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Table 3. Initial injection of an incident fast 
shock via parameter perturbation at the 
inner boundary of 25 solar radii. 



type 


slow solar wind 


location of shock nose 


if. = 0° 


width of shock front 


12° 


temporal duration 


1 hour 


shock speed 


1630 km s" 1 


ratio of total pressure 


24 



Table 4. Maximum values of parameters (z 2 G, z 2 Gr, z 2 Gt, ti) along LOSs L1-L7, used fori 
normalization in Figures [4] and [5] Each LOS is designated a time t, longitude ip, and elongation 
e in columns 2-4, and is superimposed on Figure [3] 



LOS 


time 
t 

(hour) 


longitude 

<P 
(°) 


elongation 

£ 
(°) 


intensity 
i, i R , i T 


Thomson-scattering 
geometry 
factors 
z 2 G, z 2 G R , z 2 G T 


electron 
number 
density n 
(cm -3 ) 


LI 


12 





17.5 


2.87 x 10~ 27 


2.99 x 10" 29 


104.8 


L2 


25 





40 


2.33 x 10~ 28 


6.54 x 10" 30 


35.7 


L5 


37 





90 


4.65 x 10~ 29 


2.7 x 10" 30 


17.3 


L6 


37 





155 


8 x 10~ 29 


4.92 x lO -30 


17.7 


L3 


25 


-45 


27 


4.96 x 10~ 28 


1.31 x 10" 29 


41.5 


L4 


25 


-45 


51 


1.46 x 10~ 28 


4.47 x 10" 30 


32.7 


L7 


37 


-90 


52 


1.11 x 10~ 28 


4.35 x 10" 30 


25.6 



Figure Captions 

Figure 1 (a) The Sun as a surface light source, (b) configuration for the line-of- 
sight (LOS) integral calculations, and (c, d, e) typical Thomson-scattering 
geometries. In panels (a, c, d, c), points S, Q, and O indicate the Sun 
center, a scattering point, and an observer, respectively. The solar surface 
is denoted by a grey semi-circle in panel (a). The ZSOQ is defined as 
elongation e, ZQSO as longitude <p, and ZSQO as scattering angle \- 
X < 90°, x = 90° j an d X > 90° correspond to backward, perpendicular, 
and forward Thomson-scattering, respectively. The Thomson sphere is the 
locus of points corresponding to perpendicular scattering. Point T is the 
tangent point on the solar surface from point Q, such that ZQTS = 90°. As 
viewed from point Q, the Sun has an angular half-width of ZSQT = Q, and 
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appears as a luminous disk. As shown in panel (b), for an idealized detector 
with a surface area S A and a beam size 6u>, the recorded brightness in each 
pixel of the detector is an integral through a cone determined by the size 
of the pixel's point spread function. In panel (b), a local coordinate system 
is dchned such that its z-axis is along the LOS and its xy~plane is in the 
plane of sky. All diagrams in this figure are modified from Figures 3, 6, 
and 8 of |Howard and Tappin| [20091 The nose of the interplanetary shock 
studied in this paper is initially at longitude <p = 0°. The longitude ip of 
an observer situated to the west (east) of the shock nose is defined to be 
positive (negative). Elongations e west and east of the observer's LOS are 
positive and negative, respectively. 

Figure 2 Thomson-scattering geometry for three observers at fixed longitudes 
in the ecliptic. The three observing sites, shown as white solid tiny circles, 
are at <p> = 0° (panel (a)), at p — —45° (panel (b)), and at <p = —90° 
(panel (c)), respectively. In each panel, the ecliptic cross-section of the 
corresponding Thomson-scattering sphere is depicted as a dotted circle. The 
background image shows the initial distribution of electron number density 
no- Pink and black solid lines represent the sunward and anti-sunward 
interplanetary magnetic field lines, respectively. 

Figure 3 Relative electron density enhancement (n — no) /no (a, c, e) and radial 
velocity v r (b, d, f) for an interplanetary shock propagating in the ecliptic 
at three fixed times, corresponding to 12, 25, and 37 hours. A solid arrow 
denotes the LOS from an observer located at a specific longitude tp at a 
radial distance of 1 AU. Four LOSs (LI, L2, L5, and L6) are viewed from 
<p = 0°, two LOS (L3 and L4) from <p = -45°, and one LOS (L7) from 
<p = —90°. For each of the three observers, the corresponding Thomson- 
scattering sphere is indicated as a dotted circle. Here, the elongations e 
corresponding to LOSs L1-L7 are 17.5°, 40°, 27°, 51°, 90°, 155°, and 52°, 
respectively. 

Figure 4 Intensities (i, zr, z't), Thomson-scattering geometry factors (z 2 G, 
z 2 Gr,, z 2 Gt), electron number density n, polarization p, depth z plotted 
as a function of modified scattering angle \* — 90 — X° along LI (Column 
(A)), L2 (Column (B)), L5 (Column (C)), and L6 (Column (D)). These are 
the LOSs that correspond to an observer at ip = 0° , so the shock is heading 
toward the observer, i, z'r, i-p, z 2 G, z 2 Gr, z 2 Gt, and n are all normalized 
to their respective maximum values along each LOS, as given in Table SJ 
Note that i = nZ 2 G, z'r. = tzZ 2 Gr, z't = n Z 2 Gt- For i and n, initial and 
disturbed profiles are depicted as solid and dashed lines, respectively 

Figure 5 Intensities (z, zr, z't), Thomson-scattering geometry factors (z 2 G, 
z 2 Gr, z 2 Gt), electron number density n, polarization p, depth z plotted 
as a function of scattering angle x* = 90 — x° along L3 (Column (A)), L4 
(Column (B)), and L7 (Column (C)). LOSs L3, L4, and L7 correspond to 
ip = —45°, —45°, and —90°, respectively, so the shock is propagating off 
the Sun-observer line, i, z'r,, «t, z 2 G, z 2 Gr,, z 2 Gt, and n are normalized to 
their respective maximum values, as given in Tabled 
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Figure 6 The normalization factors I*, 7^, and I T for Figures [7] and |U The 
electron number density n is assumed to be dependent on heliocentric 
distance r according to the expression n oc r~ 2 . Assuming that n oc r~ 2 and 
using the known values of the Thomson-scattering geometry factors (z 2 G, 
z 2 Gr, z 2 Gt), intensity profiles (i, ir, %t), calculated along an elongation 
of e = 30°, are presented in column (A). Integrated intensities (7*, 7^, I T ) 
as a function of elongation between 0° and 180° are presented in column 
(B). 7*, 7^, and I T are linearly scaled such that they are unity at e = 6.7°. 
In addition, the polarization distribution p is given in panel (h) . 

Figure 7 Patterns of brightness 7/7*, 7r,/7^, It /It an d polarization p, as 
viewed by observers at 1 AU and at times of t = 12, 25, and 37 hours, 
are presented in the longitude-elongation (if — e) parameter space. The 
dotted lines correspond to 7/7* = 3.68 x 10 -15 in panels (a-c), Ir/Ir = 
0.68 x 1CT 15 in panels (d-f), and 7 T /7* = 2.93 x 1CT 15 in panels (g-i). Note 
that the normalization factors, 7*, 7j^, and 7J, are shown in Figure [5] 

Figure 8 Patterns of brightness 7/7*, 7r/7^, It/I t , and polarization p, as 
viewed over the elongation range e — 6.7° - 180° by observers at 1 AU 
and at longitudes of ip = 0°, —45°, and —90°, are continuously recorded 
during the time interval of t = - 90 hours. The dotted lines correspond to 
7/7* = 3.68 x 10~ 15 in panels (a-c), 7 R /7* = 0.68 x 10" 15 in panels (d-f), 
and 7 T /7^ = 2.93 x 10~ 15 in panels (g-i). These dotted lines bound the 
brightness patterns associated with the white-light imaging of the travelling 
shock in the e — t parameter space. Note that the normalization factors, 7*, 
7^, and 7|., are shown in Figure [6l 
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(a) Solar Surface Light (b) Integrated Lines of Sight 




(c) Backward Scattering (d) Perpendicular Scattering (e) Forward Scattering 




Figure 1. 



Three observers fixed at different longitudes 



tfi = 0° <p = -45° <f - -90° 




-3-2-10 1 2 3-3-2-10 1 2 3-3-2-10 1 2 3 
r (10 1 K.) r (10" R,) r (10 ! K 3 ) 



Figure 2. 



i = 12 hours t = 25 hours t = 37 hours 




-3 -2 -1 1 2 3 -3 -2 -1 1 2 3 -3 -2 -1 1 2 3 
r (10' fl s ) r (10 ! fls) r (10 1 R s ) 



Figure 3. 
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(A) t = 12 hours (B) t = 25 hours (C) t = 37 hours 




-90 90 -90 90 -90 90 

f C) f (°) f (") 

Figure 7. 
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